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Autoionizing resonances that arise from the interaction of a bound single-excitation with the 
continuum can be accurately captured with the presently used approximations in time-dependent 
density functional theory (TDDFT), but those arising from a bound double excitation cannot. In 
the former case, we explain how an adiabatic kernel, which has no frequency-dependence, can yet 
generate the strongly frequency-dependent resonant structures in the interacting response function, 
not present in the Kohn-Sham response function. In the case of the bound double- excitation, we 
explain that a strongly frequency-dependent kernel is needed, and derive one as an a posteriori 
correction to the usual adiabatic approximations in TDDFT. Our approximation becomes exact 
for an isolated resonance in the limit of weak interaction, where one discrete state interacts with 
one continuum. We derive a "Fano TDDFT kernel" that reproduces the Fano lineshape within the 
TDDFT formalism, and also a dressed kernel, that operates on top of an adiabatic approximation. 
We illustrate our results on a simple model system. 



I. INTRODUCTION 

Since the birth of quantum mechanics, the study of 
photoionization has been important tool in characteriz- 
ing the electronic structure of materials. It is desirable 
for theoretical methods to supplement, support, inter- 
pret, and even predict the experimental spectrum. Reso- 
nance structures, arising from the interplay of bound and 
continuum excitations (where this classification refers to 
some zeroth order model) , create a fascinating panorama 
of peaks in the spectrum, whose profiles contain much 
information about the electronic states. There are three 
major issues to be surmounted in the theoretical treat- 
ment of resonances: first, resonances require an adequate 
treatment of electron correlation. We shall return to this 
point shortly, but note first that the considerable ad- 
vances in electronic structure methods and codes for exci- 
tations over the years are predominantly set up for bound 
states, not continuum states. Herein lies the second is- 
sue which is adapting the many-body methods for non- 
square-integrable scattering- type states [ll, Q . Both basis 
set issues as well as the finite matrix-based algorithms 
established in quantum chemistry codes need to be re- 
visited. When the system of interest is a molecule rather 
than an atom, a third ingredient compounds the prob- 
lem: treating the correlated continuum states in multi- 
center non-spherical potentials. A variety of theoreti- 
cal methods have been developed to treat these issues 
to a variety of extents, and we mention only a smatter- 
ing of these here. For atoms, one of the more success- 
ful approaches is multi-configuration Hartree-Fock Q; 
this accounts for electron correlation and core polariza- 
tion accurately enough to describe atomic resonances eg. 
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in halogens For molecules, Ref. [5] studied the re- 
lationship between interatomic distances and resonance 
positions using a minimal-basis static exchange method. 
Another approach utilizes R-matrix theory within multi- 
channel quantum defect theory, eg. for the calcium atom 
in Ref. Q. Reformulating the scattering problem as a 
bound-state problem in this way means that advanced 
electronic structure codes may be used. A complex- 
scaled full-configuration-interaction method was used in 
Ref. to calculate the resonances in a two-electron 
quantum dot. Methods using a complex absorbing po- 
tential in a configuration-interaction calculation, or with 
a correlated independent particle potential have been de- 
veloped; an inner-valence autoionizing resonance of the 
neon-dimer NeJ Q, the nitrogen dimer N2 and acety- 
lene C2H2 @, for example, were computed in this way. 
Several other works studied autoionizing resonances in 
acetylene, e.g. Ref. [l^l used configuration interaction 
within the multichannel Schwinger variational method, 
and Ref. [llj used Gaussian-type orbital basis sets. 

Accounting for electron correlation becomes an in- 
creasingly Herculean task for wavefunction-based meth- 
ods as the number of electrons in the system grows. 
Time-dependent density functional theory (TDDFT) 
fisl ITqI . l20j is a particularly efficient approach to the 
many-body problem, making it attractive for calculat- 
ing photo-ionization spectra of chemically interesting sys- 
tems. Recent work by Stener, Decleva, Fronzoni, and 
co-workers, Refs. [H, M, Q El, H [13 (and references 
therein) , has shown that TDDFT predicts accurate reso- 
nance parameters for a range of medium-size molecules, 
including acetylene, carbon monoxide, silicon tetrafluo- 
ride, and sulphur hexafluoride. Their earlier works used 
a one-center expansion B-spline basis, while the later 
ones utilized multi-centric B-spline basis functions, more 
suited to larger molecules. 

Due to its favorable system-size scaling, TDDFT has 



2 



become a method of choice for the calculation of bound 
spectra in quantum chemistry. For example, it allows cal- 
culations of spectra of systems as large as biomolecules 
(see eg. Refs. [H, 113, Ull), and coupled electron- ion 
dynamics on non-trivial chemical reactions (26| . One 
first computes the Kohn-Sham (KS) spectrum, which 
is the response of the non-interacting KS system: the 
single-particle excitations and oscillator strengths of the 
ground-state KS potential. The true spectrum is ob- 
tained by applying the linear response TDDFT exchange- 
correlation (xc) kernel (see also Sec. IIII|) : operating via 
a matrix equation or a Dyson-type integral equation, the 
kernel mixes the single-excitations of the KS system, and, 
were the exact kernel known, this would yield the exact 
spectrum of the true interacting system. In practise, ap- 
proximations are needed for the kernel, as well as for 
the ground-state KS potential out of which the KS spec- 
trum is calculated. Of particular note for the present 
paper, is that almost all calculations use an adiabatic 
approximation to the kernel, meaning one that has no 
frequency-dependence. 

In the TDDFT calculations of auto-ionization cited 
above, the role of the channel coupling is very clear: the 
bare KS spectra are smooth and relatively featureless, 
while after the TDDFT procedure is applied, resonances 
are generated. In one of these earliest TDDFT calcu- 
lations (in the Ne atom it was noted that while 
the resonances arising from bound single KS excitations 
whose energy lie in the continuum (eg. core to Rydberg 
excitations) are quite accurately predicted, those arising 
from bound double-excitations, are totally missed. (This 
was also noted to occur in acetylene [l^)- This was ex- 
plained in Ref. 12] as arising because the linear response 
method involves only first-order changes in the density 
(or wavefunction) , and therefore only single-excitations 
can be obtained. However, in principle, TDDFT linear 
response reproduces all excitations of the system, which 
may be linear combinations of Slater determinants with 
any number of excited electrons. The lack of resonances 
from double-excitations is not a failing of TDDFT, but 
rather is a failing of the approximation for the xc ker- 
nel that is used. For bound double-excitations, Ref. [27| 
showed that the exact kernel is necessarily strongly- 
frequency-dependent in the neighbourhood of a state of 
double-excitation character. An approximate frequency- 
dependent kernel was derived there to account for double- 
excitations; this was successfully tested on real molecules 
in Refs. [28, 29]. 

In the present paper, we investigate the form of the 
kernel that is needed in order to capture autoionizing 
resonances arising from a bound double-excitation. Our 
derivation essentially adapts Fano's 1961 analysis [lOl 
to the case when the "unperturbed" states in his con- 
figuration mixing are the relevant KS bound-state and 
the continuum its energy lies in. We begin therefore, in 
Sec. ini with a brief recapitulation of Fano's formula. In 
section Hm we discuss the implications of Fano's formula 
for the density-density response functions of TDDFT, ex- 



plaining with an illustrative sketch, that while a strong 
frequency-dependence is required in the TDDFT ker- 
nel to capture resonances arising from a bound double- 
excitation, resonances arising from bound single excita- 
tions are captured by the usual adiabatic (non-frequency- 
dependent) TDDFT kernels. Then in section [TVl we de- 
rive an approximation, in the spirit of Fano, which does 
capture the double-excitation resonance, and illustrate it 
on a simple model system. Our approximation becomes 
exact in the limit of weak interaction, for an isolated nar- 
row resonance in a single continuum. 

II. FANO RESONANCES 

The universality of the Fano profile has been noted by 
many, from lineshapes in spectra of atoms, molecules, 
solids and heterostructures, to interference in quantum 
dots and Aharanov-Bohm rings; its robustness reflected 
in the more than 3000 citations of his 1961 paper, 
Ref. [lOl- The simplest type of resonance occurs when 
a single bound-state interacts with a single continuum. 
In the absence of this interaction, the continuum is as- 
sumed to be relatively "flat", i.e. featureless. Fano per- 
formed a careful diagonalization of the Hamiltonian for a 
continuum coupled to a single bound state j^b) whose 
energy lies in the continuum [s^l- Denoting the cou- 
pling Hamiltonian by Vcp\, one defines the matrix ele- 
ment, Ve — {^E\Vcpi\^h), with being the (uncou- 
pled) continuum state. Fano derived the following for- 
mula for the matrix element of some transition operator 
T (eg. a dipole operator) between an initial (bound) state 
\i) and a state lying in the resonance region, resulting 
from the diagonalization: 

\{^E\nw (^ - c^.)2 + (r/2)2 

Atomic units are used throughout this paper. Here uj = 
E — Ei is the frequency, and 

ujr = Er-E, = E^-E,+P J Y^dE' (2) 

is the "position" of the resonance, shifting the unper- 
turbed bound-state frequency -Bb — Ei by the principal- 
value integral. The parameter 

T^2n\VEf (3) 

defines the width of the resonance, while the parameter 
q characterizes its asymmetry: 

_ (j>b|r|») + PjVE'{^E'\f\t)/{E - E')dE' 
nVE{'i>E\f\i) 

For example, q = represents a negative purely sym- 
metric Lorentzian, q —>■ oo represents a positive purely 
symmetric Lorentzian, and q = ±1 represents a purely 
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antisymmetric lineshape. The asymmetry can be inter- 
preted as a consequence of interference between the au- 
toionizing state and the continuum states [13, HH, H^l ■ In 
Refs. [13, [13, it is argued that typically q is negative. 

Although Eqs. ^ and ^ appear energy-dependent, 
q, r, and LUr are regarded as constant through the res- 
onance region. For a narrow enough resonance, this is 
a reasonable assumption; essentially the idea is that T is 
the smallest energy scale in the system. Fano also derived 
a "sum-rule" for the integrated transition probability: 

\{^^\mf = I dE{\{^E\mf-\{<i>E\mf) 

= \{^E\m\'^{q'-l)T (5) 

which expresses the unitary nature of the diagonaliza- 
tion procedure. Eq. [5] is essentially a consequence of the 
following closure relation: 

J \-^E){'^E\dE = J \^E){^E\dE+\<^^){^^\ (6) 

Fano's Eq. ([T]) tells us that the transition to the con- 
tinuum when a discrete state couples to the continuum, 
is equal to that without the coupling, multiplied by a 
(generally asymmetric) Lorentzian line-shape factor. 

Photoabsorption cross-sections measure the dipole 
transition probability, so take T to be the dipole oper- 
ator. Fits are routinely made for the Fano parameters 
q, r and LOr for a given cross-section obtained from ex- 
periment or theory, i.e. Eqs. ([2] - HI) are not typically 
used to calculate these quantities, rather, they are ex- 
tracted from experimental or theoretical data. Although 
we will only use the simplest type of resonance in the 
present paper, we do note that Fano's analysis has been 
generalized in several directions [13, 112, [13|' ^-S- when 
more continua are present. The more complicated situa- 
tions involve more fitting parameters, and in practise fits 
are made for these more general formulae, rather than 
the simplest situation discussed above. 

III. AUTOIONIZING RESONANCES WITHIN 
ADIABATIC TDDFT 

A. Photoabsorption/ionization in TDDFT 

The photoabsorption or photoionization cross-section 
measures the linear density response 6n(T,uj) of a sys- 
tem subject to an external electric field of frequency us. 
Letting a, (3 denote the three spatial directions, the cross- 
section tensor is 

CTa/3(t^) = ^y" dVdW„r^/mx(r,r',cj) (7) 

where x(r, r', t — t') = 6n{r, f)/5woxt(r', t) is the density- 
density response function of a system to an external po- 
tential ?;cxt(r, (see, eg. Refs. [13, H^)- A sum-over- 
states expression may be obtained using the standard 



linear response theory expansion for the density-density 
response function: 

_ ^ {0\h{r)\E){E\n{r')\0) {0\n{r')\E){E\h{r)\0) 
' ' ' ^ Lj-{E-Eo) + iO+ Lu + {E-Eo) + iO+ 

E 

(8) 

where \E) label the excited states, and |0) is the ground- 
state with energy i?o. n{v) = 5{y — r;) is the density 
operator. Assuming real eigenstates, the imaginary part 
is extracted as 

3x(i-,r',w) = -7r^(0|n(r)|£;)(£;|n(r')|0) X 

E 

{8{uj-uje)-5{oj + oje)) (9) 
where uje = E — Eq. Inserting this into Eq. [71 we have 

a^.iu;) = -l^^(0|f„|£;)(£;|r;3|0)x 
^ e 

{S{uj - uje) - S{uj + uje)) 
= -^^^y^^da{Lj)di3{uj) {S{uj - uje) - 5{lj + 

^ E 

where ^0.(0;) is the transition dipole moment from the 
ground-state to the excited state of energy E — uj + Eq. 

Therefore, to compute photoabsorption or photoion- 
ization, one needs an efficient way to calculate either the 
density response 6n(r, uj), or the density-density response 
function, x(r,r',Ci;) (Eq. [71) or the excited states of the 
system (Eq. [T0| . Given that the electrons in the system 
are interacting with all the others via Coulomb repulsion, 
this becomes a daunting task for correlated wavefunction 
methods for all but the smallest molecules: the numeri- 
cal effort in solving the problem scales exponentially with 
the number of electrons in the system. 

In TDDFT, instead of dealing with the correlated 
many-body wavefunction, one solves for the much sim- 
pler single-particle orbitals (r, t) that evolve in the one- 
body KS Hamiltonian. Nevertheless, from the (f>i{r,t) 
in principle, not just the exact density but all proper- 
ties of the true interacting correlated system may be ex- 
tracted [H, 1131 . In practise approximations are needed 
for the xc terms. 

Linear response TDDFT [2l|, |23| is founded on the re- 
lation between the interacting and KS density-response 
functions, both of which are functionals of the ground- 
state density of the system of interest, no(r). First, we 
note that for non-interacting systems such as KS, the 
numerator in Eq. [51 simplifies to products of occupied 
(indexed by i) and unoccupied (indexed by a) orbitals: 

. _ V $.a(r)$*JrO <^>Ur)<^^a{r') 

' ' ' ^ - (e, - £,) + zO+ UJ + (e, - e,) + iO+ 

I, a 

(11) 

where $ia(r) = 0*(r)0Q(r). The orbital energy- 
differences in the denominator — , are the KS single- 
excitation frequencies, and the orbitals live in the KS 



4 



potential Vs{r) = Wcxt(i") + Wh(i") + Wxc(i"). Here Vair) = 
/ d'^r'n(r')/|r — r'| is the Hartree potential, and Wxc(r) 
is the (ground-state) xc potential. The frequencies of 
the true system lie at the poles of the true interacting 
response function, x(r,r',cij): 

x[no]{r,r',uj) = Xs[no]ir,r' ,uj) + 

J dridr2Xs K] (r, ri , w)/hxc ["o] (^i , ^2 , w)x[no] (r2 , r', i4)12) 

where /hxc denotes the "Hartree-exchange-correlation 
kernel" , 

/Hxc["o](r,r',t^) = 1 — ^ + /xcM(r,r'w) 
|r — r' 

= X7'[noKr,r',Lu)-x-'[no]{r,r',uj) (13) 

with /xc[no](r,r',i - *') = Svxc[no]{r,t)/Sn{r' ,t). 

Although the exact xc kernel, /xc('^) is frequency- 
dependent, reflecting a dependence on the history of the 
density in the time-domain [3^, the majority of calcula- 
tions today utilize an "adiabatic" approximation, mean- 
ing one where Wxc [^] ('": t) depends only on the instan- 
taneous density /^^[no] (r, r', t — t') oc S{t — t'). When 
Fourier-transformed, this gives no structure in frequency- 
space. Instead, the adiabatic approximation is based on 
a ground-state energy functional 

\nn](r r') = ^^^-^i""^ (14) 

Despite incorrectly lacking frequency-dependence, the 
adiabatic approximation yields remarkably accurate re- 
sults for most excitations. Eq. [T2]can be transformed to 
a matrix equation indexed by KS single-excitations [2^ . 
the workhorse of bound-state TDDFT. When the true 
interacting state is composed of mixtures of KS single 
excitations, adiabatic TDDFT is expected to work rea- 
sonably well - provided the spatial functional depen- 
dence of Vxc(j) is adequately nonlocal for the prob- 
lem at hand. For states of multiple-excitation charac- 
ter however, one must go beyond the adiabatic approx- 
imation as was shown for bound states of finite sys- 
tems [13, m, [H, [13, [HI , as mentioned in the introduc- 
tion. We shall show shortly, that when such a state lies 
in the continuum, a frequency-independent kernel misses 
the resonance arising from it. 



B. Autoionizing Resonances in TDDFT 

In photo-ionization, the energy of interest lies in the 
continuum; the sum over the delta-function peaks in 
Eqs. ([9]) and (flO)) becomes an integral, and, (for positive 
to, greater than the ionization threshold), 

9x(r,r', w) = -TT (0|n(r)|£; ^ lo + Eq) {E uj+Eo\h{r')\0) 

(15) 



where the continuum states are chosen real and energy- 
normalized. The KS cross-section for continuum transi- 
tions are largely structureless, especially when only one 
continuum is relevant, the KS transition dipole is: 

ds(a;) = ((/),=,,+^|r|0,} (16) 

where (pi is the KS occupied orbital, of energy e^, out 
of which excitation occurs to the continuum orbital (j)^, 
of energy e = ei + uj. This typically gently decays as a 
function of frequency, reflecting the decay of the overlap 
between an occupied orbital and a continuum one as the 
energy of the continuum state rises. 

When one applies TDDFT to obtain the spectrum of 
the interacting system, the kernel /hxc mixes the KS 
single excitations and the spectrum distorts to varying 
degrees: less at the higher frequencies, and most signifi- 
cantly near resonances. The TDDFT kernel smears the 
oscillator strength from bound transitions whose energy 
lies in the continuum over a narrow range in the con- 
tinuum, implicitly performing the job of the Fano diago- 
nalization of Section [Til Effectively, a rather featureless 
KS continuum spectrum is turned into a dramatically 
frequency-dependent interacting spectrum, via the oper- 
ation of the xc kernel. Refs. [IM [ll, [H, [lE [H, [13 have 
demonstrated this explicitly on a wide range of interest- 
ing atoms and molecules. These works show that adia- 
batic kernels reproduce resonance features of interacting 
systems rather well through this action, when the reso- 
nances arise from a bound single excitation with energy 
lying in the continuum (e.g. of core-Rydberg nature). 

The appearance of the resonance can therefore be un- 
derstood as a mixing of the single excitations appear- 
ing in Eq. (|lip via Eq. (fT2)) . Yet, how a frequency- 
independent kernel can transform the largely frequency- 
independent KS continuum spectrum into a spectrum 
that does have such dramatic frequency-dependence, may 
strike one as incongruous. To understand this better, we 
first find an expression relating the interacting and KS 
response functions near an autoionizing resonance, fol- 
lowing in Fano's footsteps. 

We keep to the simple case of one discrete state lying 
in one continuum. Further, we consider, initially, only 
the resonant coupling; that is, we treat the KS states as 
Fano's "pre-diagonalized" states, without accounting for 
mixing amongst them. Although this could only be the 
case when there is no electron-interaction, the justifica- 
tion for this simplification is that in the vicinity of the 
resonance, this resonant coupling is certainly the dom- 
inant effect: the coupling amongst the KS states is an 
order of F less. (Later, in Sec IIVI we relax this assump- 
tion). For the matrix elements on the right-hand-side of 
we use Eq. ([T]), but taking T as the density operator n{r): 

where the ket on the right is the KS continuum excited 
state at energy E = uj + Eq- The state |0) on both sides 
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of the equation is the initial state out of which transitions 
occur, which we will take to be the ground-state. Techni- 
cally, this should be the interacting ground-state because 
Fano's analysis assumed everything except for the reso- 
nant coupling was in the zero-order states. However, if 
we approximate the bra (0| on the right of Eq. (fT7|) to 
be instead the KS ground-state, this formula then di- 
rectly relates the interacting matrix element (related to 
the oscillator strength) to the KS matrix element. This 
approximation holds if the KS ground-continuum transi- 
tion in the absence of the resonance is a good approxima- 
tion to the true interacting transition, and again, holds 
under the justification that the resonant coupling is the 
dominant effect in the resonance region. 

The quantities F, q and uJr are given by Eq.[31[31 andUl 
where the "coupling" Vi^pi is the difference in the Hamil- 
tonian of the true and KS systems: 

Kpl = Vcc -Vu- Wxc (18) 

where 14e = \ X^ij l/l'^i ^ '^jl ^he electron-electron in- 
teraction. Inserting Eq. [T7] into Eq. 1151 for uj near a 
resonance, 

^ _ (c^ - CO, + Fg(r)/2)(L^ - LOr + Fg(r0/2) 

[lj - LOrY + (r/2)'^ 
X 3xs(r,r',co) (19) 

where Xs(r,r','^) denotes the continuum contribution to 
the KS response function at frequency uj. Following di- 
rectly from Fano's analysis, Eq. [Tn]is a new relation be- 
tween the true and KS response functions in the neigh- 
borhood of any isolated narrow autoionizing resonance, 
that arises from the mixing of a single discrete state with 
a single continuum, under the assumption of weak inter- 
action. The real part of the response function can be 
obtained via the principle- value integral: 

^X(r,r-,.)^^P r ^'^f^'^'f d.' (20) 

due to known analyticity properties of x (eg- Ref. [111). 

We now address the curiousity raised earlier regard- 
ing how application of a frequency-independent adiabatic 
kernel to the flat KS spectrum can generate frequency- 
dependent resonant structure in x- A simple sketch is 
instructive to show this. Consider a resonance due to a 
bound single-excitation, at frequency Wf, > /, the ioniza- 
tion potential. Within our weak-interaction assumption 
the KS single excitation couples only to the continuum 
states, so, near Wb, Xs has the form 

Xs(c.) = Xs(^)+X^(c^) (21) 

where Xs(^) is complex, smooth and gently-decaying 
and Xs(^) is the contribution of the bound state to 
the sum over states Eq. (fTTT) . For the present illustra- 
tive purposes, we neglect the spatial-dependence. We 

take Xs(^) = ( + J j in the plots, although 



our conclusions in no way depend on this form, and 
Xs(w) — ^_J^^iQ+ ■ Here b represents the orbital prod- 
ucts appearing in the residue of Eq. (|lip . The imaginary 
part of this Xs (i^) is plotted as the solid line on the top left 
panel of Figure [TJ In particular, note that Q^Xsj which is 
directly related to the measured cross-section (Eq.[71), has 
the structureless continuum contribution plus a delta- 
peak (indicated by the arrow) at uj\y (taken to be 3, while 
b is taken to be -0.01 in these plots and / = 0.5au); other 
values yield similar plots). (Note that the delta-peak 
is not evident in the smooth KS cross-sections plotted 
in the graphs of Refs. [12, llH, El El |I3, because 
there only the KS continuum transitions in Xs (Eq.- (HH)) 
are included). Now consider inverting Eg. 1211 Interest- 
ingly, this immediately displays resonance-structures in 
both the real and imaginary part of x^^('^), as shown by 
the solid lines on the right-hand panel of Figure [TJ This 
may be simply seen mathematically: inverting Eq. [^l] 
reveals a Lorentizan denominator in both the real and 
imaginary parts. Applying the TDDFT kernel, to obtain 
X^^ — Xs^ — fnxc (Eq- CHI)), we see that an adiabatic ap- 
proximation /jfxc (Eq- (|14p ) just uniformly shifts the real 
part of Xs^i while not adding any additional frequency- 
dependence- Indeed the resonance structure of xi"^ and 
the adiabatically shifted one (x"^)"""^ resemble that of the 
Fano profile of Eq. p9)) : for a constant q — —3, the 
inverse of x computed from Eqs. (fT^H^O)) is plotted in 
the inset on the right. No frequency-dependence is re- 
quired in the kernel itself to obtain resonances in Xs- 
We now invert (x"^)"^ to obtain (x^), plotted back as 
the dashed curves in the left-hand panels of Figure [T] 
These resemble the Fano profile of Eq. for q = —3 
plotted in the insets on the left panels. The adiabatic 
shift in the inverse response functions, had the effect of 
Lorentzian smoothing the delta- function peak in Sxs, 
producing a resonance in 9x^- Therefore, making a adi- 
abatic shift in the real part of the inverse response fn, 
SRxs"^: turns, upon inversion, a delta-peak of 3xs into a 
Lorentzian resonance peak of ^x^- In this way, inver- 
sion of Xs^y shifted frequency-independently, generates 
frequency-dependent structure in x- 

For a resonance arising from a double-excitation how- 
ever, the analogous figures clearly show that frequency- 
dependence in the kernel is absolutely required. The KS 
response function of Eq.l^now consists solely of the gen- 
tly decaying first term, Xsj (left-hand panel of Figure 
i.e. there is no bound-state contribution to Xs because a 
double-excitation has zero oscillator strength [23|. The 
true response function however displays a resonance with 
an antisymmetric imaginary part (q = ±1) in order to 
preserve the oscillator strength sum-rule (see Sec. Ill| . 
On the right panel of Figure [2] we invert this smooth Xs 
to obtain the featureless x^^ shown (solid lines). The 
dotted curves represent the true x~^j computed using 
Eq. (jl9p for q ^ I. A frequency-independent kernel again 
only shifts the real part of Xs^ uniformly as shown, but in 
this case it cannot generate the resonance-structure of the 
true x^^i shown as dotted lines in the figure. For that. 
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FIG. 1: The imaginary and real parts of Xs(w) (left-hand pan- 
els, solid lines). The delta-function in the imaginary part of 
Xs corresponding to the bound-state, is shown as an arrow. 
On the right, we show Xs^ (solid) and the shift that an adia- 
batic kernel /hxc produces, yielding (x^) (dashed). When 
inverted to yield the adiabatic response kernel x'*(<^), the res- 
onance emerges, as indicated by the dashed lines on the left 
panels. The adiabatic curves should be compared with the 
insets, which come from Eq. [19] with q — —3, see text 



a frequency-dependent kernel of the form derived in the 
next section is required. The left-hand panel shows the 
real and imaginary parts of Xs and x"^, lacking resonance. 
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FIG. 2: The imaginary and real parts of Xs('^) (left-hand 
panels, solid lines). On the right panel, we show Xs^^ (solid) 
and {x^)~^ (dashed) that arises from application of an adi- 
abatic kernel /hxc b.s indicated. The dotted lines show the 
true inverse kernel x^^ (from Eq. [19] with q — 1, see text) 



IV. DRESSED TDDFT FOR RESONANCES 
ARISING FROM DOUBLE EXCITATIONS 

We are interested in the frequency-range near a res- 
onance that arises when the energy of a bound double- 
excitation lies above the single-ionization threshold. As 
argued in the previous section, the KS response function 
displays no resonance, and an adiabatic kernel cannot 
generate one. In this section we derive the form of the 
kernel that is required in order to capture this kind of res- 



onance, based on Fano's approach. We again stay within 
the assumption of weak interaction, where the dominant 
coupling between the KS states near the frequencies of 
interest, is the resonant coupling of the bound double- 
excitation with the continuum states. We again assume 
an isolated resonance: one discrete state coupling to one 
continuum. 

Eq. [THl relates the imaginary part of the true response 
function to that of the KS response function via the Fano 
lineshape. Our job is now to use Eq. [13] to find the im- 
plied structure of the xc kernel. 

A first simplification is that q{r)'^ — 1. This follows 
from the sum-rule Eq. [5] Within the assumption that 
\i) can be approximated by the KS ground-state, then 
the RHS must be zero, since T is a one-body operator 
and <i>b differs from the KS ground-state by two orbitals. 
That = 1 is consistent with the oscillator strength 
sum- rule (eg. Ref. [M]): both the KS system and the in- 
teracting system satisfy the oscillator strength sum-rule. 
As the KS spectrum does not contain the resonance aris- 
ing from a doubly-excited state, this suggests that the 
integrated area under the Fano lineshape factor should 
be zero. That is, the line-shape should be antisymmet- 
ric, so q{r) should be ±1. 

Examination of Eq [3] reveals that q must be positive 
for our case: the first term in the numerator on the 
right is zero by the above argument, while the second 
term divided by the denominator is positive (30l | , because 
{^E'\T\i) grows larger than {^E\T\i) where E - E' > 0. 

So we may conclude for the case of a double-excitation 
resonance, in the weak-interaction limit, that the imagi- 
nary part of the response function is given by: 



3x(^) = 



(c. 



(22) 



where T and uJr are given by Eqs. ^ and with Vcp\ 
given by Eq. (|18p . An immediate implication is that the 
spatial-dependence is unchanged: this is true only within 
the approximations stated above. 

Returning to our search for /hxc 7 Eg. [T3] requires that 
we invert the response function. For this, we need to 
first calculate its real part, which we may obtain from 
the principle- value integral, Eq. ([^0]) . Subtracting out 
3Xs from 3x, we write 



2r 



-p 



(a;'2-w2)((^_^^)2_^(r/2)2) 
(23) 

In the spirit of the Fano analysis, we assume that the KS 
orbitals are slowly-enough varying in frequency, that $Jxs 
can be pulled out of the principle-value integral. That is, 
that is relatively flat in the region of the resonance; 
further away from the resonance the integral vanishes as 
the lineshape factor decays rapidly. We obtain 
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-(r/2)^ + jc' + co^.)) (1 + 1 tan-i(2c../r)) - ^^co' - u:^ (r/2)^) In (ij^f^) 
((c.-c.,.)^ + (|)2) ((c. + ...)2 + (E)2) 



_i24l 



Consistent with our assumptions, we take F << (ujr — I), 
and considering frequencies w near ujr, we obtain, after 
some algebra. 



rV2 



(CO - c.,)2 + (£)2 



5Xs(to) (25) 



Putting Eqs. [55] and [21] together, we have 



X = Xs + 



T{T/2 + i{uj 



5Xs 



(26) 



The complex lineshape on the right relates the interacting 
response function to the non-interacting one with a dra- 
matic resonance structure. The cross-section obtained 
from the imaginary part, reproduces the Fano formula. 

One should not be alarmed by the poles in our approx- 
imate X in the upper half plane, given that the exact x 
should be analytic in the upper half plane, and its inverse 
X~^ analytic for 3(a;) > 0. But our approximate kernel 
holds only for real frequencies, moreover, for frequencies 
in the restricted range near the resonance. 

We now use Eq. to extract the frequency- 

dependent kernel. Subtracting the inverse of Eq. [25] 
from that of fEg. [Ojl we find the Hartree-exchange- 
correlation kernel to be 



/hxc(^^) 



Xs'M- xs(t^) + 



T{Tl2 + i{ijj -ujr)) 



SXs(w) 
(27) ' 

In Eq. [27l the spatial-dependence has been omitted for 
clarity, i.e. all arguments of /hxc and Xs are understood 
to be (r,r',Ci;). We call this the Fano TDDFT kernel, in 
that it reproduces exactly the Fano lineshape relative to 
the KS spectrum, when a bound double-excitation lies 
in the continuum. Its frequency-dependence is essential, 
as demonstrated by the sketches of the previous section. 
For illustration, we plot the real and imaginary parts of 
this in Figure [3] where we take Xs('^) ~ , + , ^ , as 
in the earlier plots, F = 0.1 and ujr — 3. 

This kernel accounts only for the coupling of the KS 
bound state with the continuum, which is the dominant 
effect near the resonance. It becomes exact in the limit of 
weak interaction, for an isolated and narrow resonance, 
where F is the smallest energy scale of the system. 

We may also include this kernel on top of an adia- 
batic kernel, in order to account also for mixing of non- 
resonant single excitations that adiabatic TDDFT may 
capture well. For the case that the resonance arising from 
the double-excitation is well-isolated from any other res- 
onances in the system, modifying Eq. 1261 we assert the 




FIG. 3: Sketch of the Fano TDDFT kernel, Eq.[27l that repro- 
duces Fano's lineshape for the autoionizing resonance arising 
from a double-excitation. 



dressed response-function approximation: 

. F(F/2 + .(^-^,.)) ^ . 

(c.-C..)2 + (£)2 



(28) 



where x^{^) is the interacting response function com- 
puted using an adiabatic xc kernel. At frequencies mov- 
ing away from the resonance arising from the double- 
excitation, our dressed response function reduces to the 
adiabatic response function, x'^i^): thus our dressed 
Response function leaves untampered the usual reason- 
ably accurate response that adiabatic TDDFT gives 
for single-excitation states. The corresponding dressed 
cross-section, obtained from its imaginary part, yields 



From Eq. [^5] we derive the dressed kernel: 



/hxc(^) = f^. 



(x") 



A\-l 



1 + ix^r'^x 



AT{T/2 + i{uj-Ur)) 
(C. - C.,)2 + (£)2 



- 1 30) 



In practise, computing cross-sections using our dressed 
TDDFT would proceed very simply: First an adiabatic 

calculation would be run, as in Refs. [13, m, H m, [n, 

\vi\ . Then where a bound-double excitation is known to 
lie (for example, by summing KS orbital frequencies), uti- 
lize the relevant KS orbitals in Eqs. ([3]) and ([2]), to find 
the width F and shift of the resonance position F{E) 
that appears in w^. Then modify the cross-section com- 
puted using an adiabatic approximation by the lineshape 
(Eq. I29p . That is, in practise, if interested in comput- 
ing the cross-section, we would not need to utilize Eq. [30] 
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directly, instead we would use Eq. (55] Instead, Eq. [3D] 
and Eq. [27| are of fundamental interest here: it is the xc 
kernel for which approximations are needed in TDDFT, 
and these equations reveal the form it requires, in order 
to reproduce the Fano resonance arising from a double- 
excitation. 

A. Model example 

We illustrate our results on a simple model involving 
two electrons in one-dimension living in an external po- 
tential of the form: 

Ve.t (x) = ~ l3J{l-t&nh\x))^ (31) 

cosh (ax) * 

We choose values of parameters Uo,a, /3 such that there 
are (at least) two bound single-particle states in the non- 
interacting problem. A sketch is shown in Figure [H A 
double-excitation to the first excited orbital is shown on 
the left; this has energy 2e where e is the energy difference 
between the single-particle orbitals. This energy 2e ex- 
ceeds the single-ionization threshold for this system, and 
lies in the continuum; therefore the state on the left is 
degenerate with a single-excitation to the continuum, in- 
dicated on the right. When electron-interaction is turned 
on, an isolated resonance is created, of the type to which 
our formula and analysis applies. 













2e 






Vl 


e/ 



FIG. 4: Sketch of model potential, indicating the two excited 
non-interacting states, one a double-excitation that is bound 
(on the left) and the other a single- excitation that is unbound 
(on the right). Allowing the electrons to interact couples these 
degenerate states, creating a resonance. 

We first find the KS potential and spectrum. We 
choose a weak delta-function interaction: 

= XSix ~x') , A < 1. (32) 

For small enough A, exchange dominates over correla- 
tion. We use the exact- exchange approximation for two 
electrons: 

v^c ^ ~Vh/2 ^ -Xn{x)/2 (33) 
The KS potential follows as: 

Vs {x) = ^;oxt +VH + Vy:c = 757 — - (34) 

cosh (ax) 



if we take /3 = A in Eq. [311 The exact one-electron eigen- 
states and energies of Vs (an "Eckart well" ) can be found 
in many quantum mechanics textbooks. The parameters 
Uo, a are chosen such that there are at least two bound 
one-particle states: we chose Uq — 1.875, and a = 1, 
which places the non-interacting bound orbital energies 
at eo = —1.125 and ei = —0.125. From the one-electron 
orbitals we calculate the KS response function Xs and 
dipole moment, shown as the dashed line in Figure [5l 
after placing two electrons in the lowest orbital. As ex- 
pected, the dipole moment to the continuum states is 
smooth and gently decaying. 




1.6 1.8 2 2.2 2.4 

CO 

FIG. 5: The square of the dipole moment: Dashed line is KS, 
dotted line is the adiabatic TDDFT (exact-exchange), solid 
line is the resonance with the Fano lineshape as would be 
obtained with the Fano TDDFT kernel Eq. [27l and the (x) 
are our dressed kernel used in conjunction with the adiabatic 
approximation, Eq. [30l 



We now consider the TDDFT spectra, and will com- 
pute everything only to first-order in the interaction 
strength, A, whose numerical value we take as 0.2 in 
the calculations. We first apply an adiabatic kernel, and 
again choose exact-exchange for this: 

/4(a;, x') = f^ix, x') = -X5{x - x')/2 (35) 

For the cases where the states are mixtures of sin- 
gle excitations, exchange effects dominate in the weak- 
interaction limit, being of 0(A) while correlation is 
O(A^); so for small values of A, is expected to be quite 
accurate for energies and oscillator strengths of states of 
single-excitation character. Then 

X'^(r,r',w) = xs(i-,i-',w)+^ j Xs(r,ri,w)x(ri,r',w)dVi 

(36) 

where, on the right, we may replace x with Xs to get 
the 0(A) adiabatic TDDFT spectrum. In one-dimension, 
from Eqs. [7| and [TUl 

\d{uj)\^ — / x'^x{Xtx' ,uj)x' dxdx' (37) 
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The KS version of Eq. [15] simplifies to 

'^Xs{x,x',Uj) = -TT(j)o{x)(j)o{x')(j)e=uj+eo{x)(t)e=u+eo{x') 

(38) 

Using Eq. [38] and its principal- value integral to get the 
real part, Eq. [37] finally gives 

1^-^(^)12 = \ds{Lu)\^ + 2Xds{uj) X 

/ dxi4{xi)Mxi)P I dsie')cp,,{xi) , „ rf<39) 

This adiabatic dipole moment is the dotted line in Fig [5] 
it smoothly shifts the KS spectrum, redistributing oscil- 
lator strength from the lower to higher frequencies of the 
KS spectrum, but is a small and smooth correction. 

Applying now the frequency-dependent kernel Eq. [27] 
reveals the resonance (solid line) , while the dressed kernel 
Eq. [30] (x points) incorporates corrections from the adia- 
batic approximation as well as capturing the resonance. 
As mentioned in the previous section, in practise, we uti- 
lize Eq. 1291 evaluating F and tOr using the KS orbitals. 

V. CONCLUSIONS AND OUTLOOK 

We have discussed the implications of Fano's resonance 
formula for the xc kernel of TDDFT. In the case of a 
narrow isolated resonance involving one discrete state 
and one continuum, and weak interaction, we derived 
the Fano- equivalent formula for the imaginary part of 
the TDDFT response function (Eq. [TU). We illustrated 
how a frequency-independent kernel applied to a largely 



frequency-independent KS response function, actually 
yields dramatic frequency-dependent resonant structures 
for the case of a bound single-excitation, but yields no 
structure in the case of a bound double-excitation. 

Within our assumptions, we derived the exact form of 
the frequency-dependent kernel that is needed for the lat- 
ter case: we call this the Fano TDDFT kernel (Eq. [27]) . 
in the sense that it exactly reproduces Fano's line- 
shapes, when applied to the KS density-response func- 
tion. We then asserted a dressed frequency-dependent 
kernel (Eq. [30| , that accounts for the Fano effect on top 
of an adiabatic approximation. The form of these kernels 
is of fundamental interest for TDDFT. In practise, we 
propose one computes the resonance width F and posi- 
tion LUr using the appropriate KS orbitals; compose from 
them the Fano lineshape (Eq. [29)l and thereby modify 
the adiabatic spectrum. How this will work in practise 
for molecules of interest, remains to be tested. 

In this paper we considered the simplest case of a res- 
onance arising from a bound double excitation: that is, 
when the resonance is isolated from all others and in- 
volves only one continuum and one discrete state. We 
have not considered the interaction of resonances [39| 
that arises in most systems, nor the computation of in- 
dividual branching ratios when several channels are in- 
volved jsj. The present work is only a first step in un- 
covering how the exact xc kernel captures resonances in 
the general case, and how to approximate it in practise. 
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